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Abstract 


An ever-increasing concern on the stratosphere ozone depletion has led to a 
worldwide ban on ChloroFluoroCarbons (CFCs) and prompted rigorous search for 
alternative refrigerants with zero Ozone-Depletion Potential (ODP) and lower Global 
Warming Potential (GWP). The prediction of phase equilibrium behavior of complex 
mixtures over broad ranges of temperature and pressure is an important problem in 
chemical process design. The objective of the present work is to predict and compare 
the VLE data for refrigerant mixtures using the following three different models (1) 
UNIFAC Method. (2) modified UNIFAC Method and (3) Peng-Robinson model with 
van der Waal one fluid-mixing rule. 

A total of thirty-one binary mixtures of different classes, have been studied using 
the above three models and the predicted results are compared with experimental data. 
It was found that the UNIFAC model is able to predict the VLE data much better than 


the other two models. 
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CHAPTER 1 


INTRODUCTION 

The estimation of phase equilibrium behavior of complex mixtures over broad ranges 
of temperature and pressure is an important problem in chemical process design as most 
of the industrial chemical processes deal with operations such as. extraction, adsorption, 
distillation, leaching, and absorption, which involve phase contacting. In particular 
distillation is the most important separation and purification process used in chemical 
industry, where vapor and liquid phases co-exist. This requires the process designer to 
have a reliable Vapor-Liquid Equilibria (VLE) data, for design of distillation columns. 
Though for a large number of binary and for a few multi-component systems the 
experimental VLE data can be found in literature, the experimental data is not available 
at all desired conditions. This necessitates the designer either to experimentally determine 
or estimate the required VLE data. But if the process of interest is to be performed over a 
wide range of pressures and temperatures, the number of experimental points required 
becomes large. Hence the available experimental data can be used to develop a model, 
which in turn can be used to predict the data at the required operating conditions. 

In accord with an increasing concern on the stratospheric ozone depletion, the 
timetable toward a complete ban on the production and consumption of 
ChloroFluorCarbons (CFCs) and HydroChloroFluoroCarbons (HCFCs) has been 
accelerated to phase them out. The study of the 'ozone hole' above the Antarctic 
established the relationship between the chlorine and ozone content in the stratospheric 
zone and dictated that the chlorine oxide (CIO') radical is mainly responsible for the 
ozone destruction. This radical is formed by chlorine atoms released by 
chlorofluorocarbons (CFCs) and hvdrochlorofluorocarbons (HCFCs). used mainl y in the 
refrigeration industry. The environmental regulations that have been imposed on all 
countries in order to tackle the problem of ozone depletion, forced the refrigeration 
industry to totally phase out CFCs in 1995. and HCFCs in the following years (2015) for 
western Europe) (Clodic and sauer. 1994). The recommended worldwide ban on many 
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low-molar mass chlorofluorocarbons has prompted rigorous search for alternatives with 
zero Ozone-Depletion Potential (ODP) and lower Global Warming Potential (GWP). 

The chlorofluorocarbons were primarily used as refrigerants, but are also important 
as blowing agents in polymer foam manufacture and cleaning solvents in electronic 
circuit manufacture. In many cases, such as refrigeration applications (which is most 
important area of use of CFC). it is necessary to develop alternative fluids with 
thermodynamics properties, which are similar to CFCs they are to replace so that existing 
equipment can be used with a minimum of modification. For this purpose first 
hydrocarbons (HFCs) and hvdrochlorofluorocarbons (HCFCs) and their mixtures were 
suggested, though now fluorinated alcohols and their mixtures HFCs and HCFCs and 
other fluids are also being considered (Sauermann et al.. 1993: Laugier et al.. 1994: 
Nowaczvk and Steimle. 1992). The use of mixtures can result in energetic improvements. 

extensions of the applications, such as flammability, insolubilitv in oil or green house 
effect implications. 

The available models for the prediction of VLE data can be broadly classified as 
activity coefficient models and Equation of State (EOS) models. In activity coefficient 
model the liquid phase non-ideality is described in terms of an activity coefficient (yfl for 
component / of the system. The vapor phase non-ideality is expressed in terms of fugacity 
coefficient (fc) for component / of the system. The fugacity coefficient can be obtained by 
making use of an appropriate equation of state. Some of the models that belong to this 
method are Margules. Van laar. UNIFAC etc.. (Rao Y. V. C. 1997). In EOS approach 
both the liquid and vapor phase non-idealities are explained by the fugacity coefficient 
where, it is indicated by (/) for the liquid phase and ff.) for the vapor phase. Some of 
the equations of state such as Peng-Robinson (PR). Soave-Redlich-Kwong (SRK). Peng- 

Robinson-Strvjek-Vera (PRSV). Benedict-Webb-Rubin (BWR) are used in this method 
(Marketal. 1991). 

The activity coefficient models give satisfactory results at low to moderate pressures. 
However, their applications at critical and supercritical pressures yield inaccurate VLE 
data. EOS approach is a powerful tool for addressing this problem. The major advantage 
of the EOS approach is its applicability over a wide range of temperatures and pressures 
including critical and supercritical, where the activity coefficient models fail. In process 



simulation calculations, the equations of state are widely used due to their algebraic 
simplicity and their accuracy. 

There are two requirements for the equations of state to be successful. (1) The EOS 
must predict the saturation pressure of pure substances accurately. (2) The mixing rules 
must be available which correctly extend these equations to multicomponent mixtures. 
Different equations of state predict the thermodynamic properties with different degrees 
of accuracy and they require different types of data for proper evaluation. Hence, a 
reliable EOS for practical applications must be physically sound, reasonably simple and 
present no convergence problems. 

Cubic equations of state, which are in wide use. model the repulsive forces 
incorrectly resulting in inaccuracies in representing the thermodynamic properties, 
especially in the supercritical fluid and condensed phase regions. Van der Waals 
proposed the first cubic equation of state in 1873. The well-documented van der Waals 
equation is the most influential EOS that incorporated both the repulsion and attraction of 
intermolecular forces in the expression. Numerous EOS have since been proposed within 
the framework of this concept (Anderko (1990)). Redlich and Kwong (1949) modified 
the attractive term of the expression of the van der Waals EOS. but retained its repulsive 
term, to satisfy the boundary conditions in the low and high-densitv limits. The Redlich- 
Kwong (RK) EOS is among the most fruitful of the van der Waals family for routine 
engineering calculations. Various modifications of the RK equation to improve its 
accuracy have also been suggested. These are represented by the Soave (1972) and the 
Peng and Robinson (1976) equations. 

The extension of equations of state to describe phase behavior for a broad range .of 
multicomponent mixtures had been more difficult due to limited applicability of the van 
der Waals (vdW) one-fluid mixing rules that are commonly used for relatively simple 
mixtures. The vdW mixing rules are incapable of representing the highly non-ideal 
mixture behavior of polar or associating fluids. Consequently, much effort in recent years 
has been devoted toward developing alternative mixin g rules to extend the cubic 
equations of state to complex mixtures [Orbey and sandier (1996); Heidemann (1996)]. 
V/ith the advent of multiparameter mixing rules, equations of state are now being used 



for the phase equilibrium calculations of complex mixtures that were traditionally 
described with activity coefficient models. 

The present work deals with prediction of VLE data for refrigerant mixtures using 
UNIFAC method. Modified UNIFAC and Peng-Robinson with Van der Waals one fluid 
mixing rule. The cubic equation has three roots. The identifications of these roots are 
very easy as. higher root correspondent to the vapor, low value corresponding to the 
liquid value and intermediate root has no significance. This makes the prediction of 
volume at different temperatures and pressures possible, by solving the compressibility 
equation of state in terms of the compressibility factor (Z) or volume. 

VLE calculations are done using UNIFAC. Modified UNIFAC and PR Equation for 
binary mixtures of refrigerants. Van der Waals one fluid mixing rule is used with EOS for 
the VLE calculations. The results are then compared with the available experimental data. 

Chapter#2 lays emphasis on the important models that are used for VLE calculations 
in the two available approaches. Chapter 3 focuses on literature survey. In the next 
chapter development of the models with mixing rule and calculation procedure is 
discussed. Fifth chapter presents the results and a discussion of the results, followed by- 
conclusions. 
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CHAPTER 2 


THEMODYNAMIC MODELS 

2.1 Equilibrium Criterion 

For a closed system containing vapor and liquid phases, with each phase consisting 
of c components, in a state of equilibrium at a constant temperature ( T) and pressure (P). 
the criterion for equilibrium is given by (Smith et al. 1996). 

H\ = P, (/' = 1.2. c) (2.1) 

where, u, and u. are chemical potentials of component / in liquid and vapor phases. 

respectively. The chemical potential u , . is given 

fi t = T, (T) + RT\n f (2.2) 

where r,(D. an integration constant, is a function of temperature only and /’ is the 

fugacity of component /' in the mixture. Using Eq. (2.2) in Eq. (2.1). the general criterion 
for vapor-liquid equilibrium reduces to 

/;=/: (/= i.2. o (is) 

The Eq. (2.3) can be rewritten as (Rao Y. V. C. 1997). 

r, x ,.C = fi y,P (i = 1-2 c) (2.4) 

where = activity coefficient of component / in the solution. 

.V, . y, = mole fractions of component / in liquid and vapor phases, respectively. 
f° = fugacily of component / in the standard state. 

(j)] - fugacity coefficient of component / in the vapor phase. 

P = Pressure at which the system is held. 

The standard state fugacity ( f" ) is the fugacily of component / at the mixture 

temperature and at some arbitrarily chosen pressure and composition. For most of the 
liquid solutions . the standard state for component /. is chosen as the liquid phase of the 



pure component / at the system temperature and pressure. With this choice of the 
standard state, f" is given by. 


f, 


0 


P, '(/>,' exp 



R T 



(2.5) 


where. ft is the fugacity coefficient of component / at saturation pressure Pf for the 

same component i at temperamre T. v' is the molar volume of liquid for component i. 

Substitution of Eq. (2.5) in Eq. (2.4) results in the basic equation for vapor-liquid 
equilibrium, which is given by 


y,x,Pfft exp 


v !(p-/r) 


RT 


= ft' y,P 


(( = 1.2. c) (2.6) 


At low pressures, the vapor phase can be assumed to behave like an ideal gas. Hence 


ft =1 and ft=\. Moreover, the Pointing correction factor (exponential term in Eq. (2.6)) 


is approximately equal to unity. At moderate pressures ft and ft are approximately 


equal and hence it is reasonable to assume that ( ft I ft ) = 1 . Thus at low to moderate 
pressures Eq. (2.6) reduces to 


r,x,P;' =y,p 

(1-2 

c) 

(2.7) 

y,P 

or v = 

' ' s 

(1-2 

c) 

(2.8) 


*,p; 

However, at high pressures these assumptions are not valid. Eq (2.3) can also be written 
as 



ft, X,P = ft y,P 


(since cp, = ( fjyft) ) 


or 

ft, x, = ft V, 

( 1 -2 ■ • 

c) 

(2.9) 

or 

V ft 

K t = — = ^~ 

x, 

(1.2- 

c) 

(2.10) 


where K- is the K factor or equilibrium constant for component i. 
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The methods available for the prediction of VLE data can be broadly classified as: 

1 . Activity coefficient Approach 

2. Equation of State (EOS) Approach 

The main difference between the two methods lies in the reference state. For the 
EOS method, the reference state is the ideal gas condition and the fugacity coefficient is a 
measure of the deviation from this state for both the vapor and liquid phases. For activity 
coefficient methods, the reference state for the liquid is the pure liquid at the same 
temperature and pressure of the mixture. The choice of the thermodynamic model to be 
used depends on 

• Nature of the components 

• The pressure and temperature conditions 

NATURE OF THE COMPONENTS 

The nature of the component decides the extent of non-ideality. Non-ideality is due 
to the interaction between the molecules, which depends on intermolecular forces. This 
interaction is more evident in the liquid phase where the molecules are packed closer. 
The choice between the EOS method and Activity coefficient method can be made based 
on the polar nature of the components. 

Hydrocarbons like ethane, propane, pentane are non-polar and can be represented 
adequately by EOS methods. Mixtures of polar compounds may be represented using 
Activity coefficient methods. Mixtures of polar and non-polar compounds can give rise 
to two immiscible liquid phases and hence a model that is capable of predicting a phase 
split must be considered. This would be an Activity Coefficient method like NRTL or 
UNIQUAC. Wilson or Van Laar models are incapable of predicting splitting of the liquid 
phase. 

TEMPERATUE AND PRESSURE 

At high pressures and at temperatures above critical, the activity coefficient models 
fail. EOS models are valid over a wide range of temperature and pressure. Since 
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consistency is inherent in the EOS approach, these are applicable even at the critical 
conditions. 


2.2 Activity Coefficient Approach 

The basis of activity coefficient approach is the activity, which is defined as 

_ f,(T.P .: Q 
' r°) 


( 2 . 11 ) 


Here a, is the activity of component i and f° is the standard state fugacity (at some 
arbitrarily defined pressure P° and some arbitrarily defined composition x,° ). The 
activity coefficient/, is defined as 

.. _ a i _ f, 


' i rO 

xj, 

The partial molar Gibbs free energy is given by 


( 2 . 12 ) 


M, -M? = § i.rcal - g.Micol ' =RT In - f ' 


f. 


ideal 


From Raoults law one can obtain 


f = v f° 

J i. ideal id i 

= RT1n^£- = RTlnr, 


r U 

x ,f, 

Using the definition of partial molar excess Gibbs free energy, one can obtain 


( 2 . 13 ) 

( 2 . 14 ) 

( 2 . 15 ) 


In / = — 
RT 


djngf] 


cn, 


( 2 . 16 ) 


T.P.n , 


where. g E = molar excess Gibbs energy. The relation between g E and the activity 
coefficients can also be expressed as 


g E =Tx,g; =RTTx, In/, 


i=i 


j=i 


( 2 . 17 ) 


The Gibbs-Duhem relation for a binarv mixture is given bv 


h.p 


ain/, ^ 


a In/, ^ 

ax, J 

— 

T.P 

l CX': J 


( 2 . 18 ) 
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Assume that g E can be expressed as a function of temperature, pressure and composition. 
That is. 

g i: =f(T.P.x,) (2.19) 

The Activity coefficient can then be derived from Eq. (2.16) and substituted into the 
following equilibrium constraint. 

y.Ptf =r l (T.P.x,)x,f l ° (2.20) 

The system of equations defined by Eq. (2.20) can be solved for the equilibrium 
temperature, pressure, and concentration. Therefore, the basis of activity coefficient 
approach is to find an accurate form of Eq. (2.19). 

The commonly used Activity coefficient models are 
r Margules Equation 
r- Van Laar Equation 
r' Wilson Equation 
r- NRTL Equation 
r- UNIQUAC Equation 
r UNIFAC Equation 
r- MODIFIED UNIFAC Equation 


2.2.1 Margules Equation 


The Margules equations, can be deduced from the expressions for excess Gibbs free 
energy given by Wohl's expansion ( Mark et al. 1991) 


g" 


RT {x ] q [ + x 2 q 2 ) 


— 






• J j | | “I” 2) Cl | -> j „ 


( 2 . 21 ) 


where 

z, = — = the volume fraction of component i (2.22) 

x,q x +x 2 q z 

q, = measure of volume of the component /. a/\ any am ... are empirical constants 
which signify two body, three body interactions. 
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Margules assumed the size of molecules to be equal (qj = q : ). After some re- 
arrangements Eq. (2.2 1 ) reduces to the Margules three suffix equation given by 


RT 


= (.4 i: x, + .4 :i x, )x, x. 


(2-23) 


where A,-_ and A : , are adjustable parameters. 

The activity coefficients y, and y : are given by 

lny, = [a u + 2(.4,, - .4,, )x, ]x; (2.24) 

lny, = [„4 :i + 2{A x: - .4,, )x, ]x," (2.25) 

The Margules three suffix equations are used for symmetrical systems ( q / = q : ) i.e. 
systems where all the components have almost equal size. 


2.2.2 Van Laar Equation 


Margules equation assumes equal size for molecules, which is not true. If the size 
difference of the molecules (q x * q : ) is incorporated into Wold's expansion. (Mark et al. 
1991) it gives 

g" 2 a^x.x.q.q. 

a__ _ — 1 = — ! ALL (2.26) 
RT x x q x + x : q : 


Simplifying Eq. (2.26) leads to van Laar equations, which are given by 

g E = 

RT x,.4 i; +x,.4 ;i 


(2.27) 


lny, = .4,, 



X] - 4|2 + X2^42| 


(2.28) 


lny. 




.4, : x, 


x,.4, ; + x ; .4 ;i 


where A/: and A;/ are adjustable parameters. 


(2.29) 


2.2.3 Wilson Equation 

For mixtures in which the components differ from each other in molecular size and 
interaction between the unlike and like molecules are different. Wilson (1964) proposed 
the following expression for the excess Gibbs free energy of a binary solution. 


i r\ 



"§- = -*i ln(x, + A, : .v ; )- .v, ln(x, + A,,x, ) 
K1 


(2.30) 


where A v _ and A 2] are two adjustable parameters, which are related to pure component 
molar volumes and characteristic energy differences. 

A , 2 and A 2 , are given by 


v. 


A,, = -^exp 


A,, = —exp 

v. 


A: A i 


RT 


a,, - /„, 


RT 


v, —molar volume of pure component /. 

/tij=energy of interaction between molecules of component / and /'. 

The activity coefficients can be derived from Equation (2.30) and are given by 

In v, = -ln(x, + A i: x, )+ x, 
lnv, =-ln(x ; + A :i x, )-x 



(2.31) 

(2.32) 


(2.33) 


(2.34) 


2.2.4 Non-Random Two Liquid (NRTL) Equation 

The Non-Random Two Liquid (NRTL) equation proposed by Renon and Prausnitz 
(1968) is applicable to partially miscible as well as completely miscible systems. 

The NRTL equation for excess Gibbs free energy is given by 


g_ 

RT 


= X \ X 2 




TA 2 


x, +G,,x ; x, + G,,x, 


(2.35) 


_ <m: ^ :: 

(2.36) 

RT 

_ _ g:\~gn 

(2.37) 

•' RT 

G r _ = exp(— a,,r l2 ) 

(2.38) 

G>, = exp(-or :i r :i ) 

(2.39) 


1 1 


where 



gij =energy of interaction between molecules of component i and j. 

The parameter a, : is introduced to take into account the non-randomness of the mixture. 
If a, : =0. the mixture is completely random and Eq (2.35) reduces to Margules equation. 
The activity coefficients can be derived from Eq (2.35) and are given by 
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(2.40) 


(2.41) 


The NRTL equation does not provide any additional advantage over Margules and Van 
Laar equations for moderately non-ideal solutions, but it represents the excess Gibbs free 
energy of strong non-ideal and partially miscible systems quite satisfactorily. 


2.2.5 UNIversal QUAsi Chemical (UNIQUAC) Equation 


To express the excess Gibbs free energy of a binary mixture. Abrams and Prausnitz 
(1975) developed the UNIversal QUAsi Chemical (UNIQUAC) model. The UNIQUAC 
equation for (g RT) contains two parts - a combinatorial part and a residual part. The 
combinatorial part takes into account the composition, size and shape of the constituent 
molecules and contains pure component properties only. The residual part takes into 
account the intermolecular forces and contains two adjustable parameters. The 
UNIQUAC equation is given by 

g h g !: g l: 

— — = — — ( combinatorial ) + — — (residual) (2.42) 

RT RT RT 


where 


( combinatorial ) = x. 


01 , „ <t>2 


In— + x, ln- 
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+ — 
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6> 9? 

Q\ x Y ln T" + <7;X:ln-r 
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(2.43) 


g_ 

RT 


(residual) = -^,x. In ( (9, + 9 Z r, , ) - q 2 x. In ( 9. + 6^ r, , ) 


^,= segment fraction of component i = x,r , /X x , r , 

v ! j ' J 


(2.44) 

(2.45) 



f \ 

#,=area fraction of component / = x, q,/Y. x ,q ; 

V / i J 

r, = volume parameter of component / 
q, = surface area parameter of component / 


r (l = ex P ( - 
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= exp 


^ a ^ 
// 


=adjustable parameter 


(2.46) 


(2.47) 


z/„= average interaction energy for the interaction of molecules of component i with the 
molecules of component /': and 
r = coordination number which is usually taken as JO. 

From Eqs. 2.42-2.47 one can obtain the activity coefficients as 

L R 

lny, = \ny,‘ (combinatorial) + lny i ( residual ) (2.48) 


lny: = \n — + — q, In - + l,- — Tx l I, 
x, 2 <p, x, ■■ 


In;/,' =7, 


l-ln| y.0,r„ 


-I 


9 j.. 

' 1/ 




(2.49) 


(2.50) 


where 


-1) 


(2.51) 


The strucmral parameters r, and q, are calculated as the sum of product of the group 
volume and parameters R k and Q k 



k 


(2.52) 


<7, 




(2.53) 


where v k '" is the number of groups of type k in a molecule of component i. The 
UNIQUAC equation is applicable to a wide variety of liquid solutions c.ommonly 
encountered by chemical engineers. 
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2.2. 7 Advantages and Disadvantages of Activity Coefficient Approach 

A simplifying assumption that is made in the development of activity coefficient 
models is that g E is independent of pressure. In view of this assumption, the activity 
coefficient approach can be applied to systems at low to moderate pressures, far removed 
from the critical. For processes that operate near or above the critical, the activity 
coefficient approach cannot be used, which is by far the greatest deficiency of this 
approach. 

The activity coefficient approach should be used to treat low pressure VLE data. The 
activity coefficient approach can be successfully applied to various types of systems: 
mixtures of alcohols with hydrocarbons, non-polar hydrocarbon mixtures and mixtures of 
alcohols with other polar compounds. In addition, systems with limited miscibility have 
been successfully treated with activity coefficient approach. 

2.3 Equation of State Approach 

In the equation of state method, a single equation of state is used to represent all fluid 
phases. From the thermodynamic point of view, this is a potentially more powerful 
approach as it provides a uniform representation of the thermodynamic properties in thfe 
two-phase and one-phase regions with a single equation. It is not only applicable to wide 
pressure and temperature ranges including critical and supercritical conditions but also 
makes it possible to calculate phase equilibria data. According to the review by 
Kolasinska (1986). depending on the form of equation of state variables, the equation of 
state may be classified into the following types. 

1 . Non analytical 

2. Analytical 

a) inspired on a Virial expansion 

b) inspired on the two-term van der Waals form 
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Equations of state inspired on the two-term van der Waals form were found to be 
most attractive due to their simplicity, low computational costs and reliability. Non- 
analvtical types of equations are not presented here. A brief account of the analytical type 
of EOS is presented below. The simplest EOS is the ideal gas equation (PV = nRT) which 
considers that there are no interactions between molecules. 


2.3.1 Virial EOS 


In real gases molecular interactions alway s exist. This is explained by the Virial 
EOS where Rl'is expressed as a power series in P (Smith et al. 1996). 

PV = a + bP + cP~+ (2.54) 

PV = a(\ + B'P + C'P' +.... ) (2.55) 

Where, a. B'.C' are constants at a given temperature for a given species. When a = 
RT. Eq. (2.54) reduces to 

Z = ^L^\ + B'P + C’P : + (2.56) 

RT 


or alternatively 


Z = 


, + « + 

V 


c_ 

V- 


D 

H - + .... 

r-' 


(2.57) 


Where. Eq. (2.56) or (2.57) is called Virial equation of state. Z is the compressibility 
factor and B. C. D etc., (or B'.C'. D'etc .. ) are called virial coefficients. The virial EOS is 
not used in phase-equilibrium studies as it can describe only vapor phase behavior. 


Equations inspired on the two term van der Waals form 

In general equations based on the van der Waals model are in a form where total 
pressure is the sum of an attractive and a repulsive pressure term (Kolasinska 1986). 


P Pi repulsive I P laitracnvei (2.58) 

This form is proposed to take into account the attractive and repulsive forces existing 
in the real fluids. Imagine two molecules as hard spheres. When the molecules are at a 
large distance from each other they experience small or no attractive forces. As the 



distance between the molecules is decreased, the attractive forces get larger, until the 
hard spheres come close to each other. Once the molecules come close to each other, 
repulsive forces start dominating until the molecules can get close no further. This 
volume is unavailable for the molecular motion. Van der Waals took the unavailable 
volume into account by introducing a parameter b called the 'excluded volume'. 


2.3.2 van der Waals Equation of State (vdW) 


When expressed in the form of Eq. (2.58) the van der Waals equation is given by. 

RT a_ 

V" 


v-b 
RT 


' repulsive 


(2.59) 

(2.60) 



( 2 . 61 ) 


The values of a and b can be found by regressing vapor pressure data and choosing a 
and b that best fit that data or alternatively critical point conditions are used. At the 
critical point, the first and second derivatives of pressure with respect to volume at 
constant temperature are zero. When these conditions are applied to Eq. (2.59). the values 
of a and b turn out to be 


21R Z T Z 

a = — 

64 P 

, RT 
8P 


(2,62) 


(2.63) 


2.3.3 Redlich-Kwong Equation of State 

An important modification of the van der Waals equation of state was 
made by Redlich and Kowng (1949). who introduced a temperature dependence and 
slightly different volume dependence in the attractive term. The Redlich-Kwong 
Equation of State (RK EOS) is given by 
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a 


(2.64) 


P = 


RT 

v-b 


T -v(v + b) 


where 

■i t ^ 

a = 0.42748—^ — (2.65) 

RT 

b = 0.08664 — - (2.66) 

P 

These modifications improved the accuracy of the equation. As before, the values of a 
and b can be found by fitting vapor pressure data or by using the conditions at the critical 
point. The RK EOS has limited accuracy, and is generally successful only for nearly 
ideal systems. To overcome these deficiencies, many modifications have been proposed 
over the years. Two of these modifications have achieved wide acceptance - the Soave 
and the Peng-Robinson equations of state. 


2.3.4 Soave-Redlich-Kwong Equation of State 


Soave (1972) proposed a modification to the RK Equation incorporating the 
temperature dependency to the EOS parameter a as 

a(T) = a c a(T ) (2.67) 

where 


S = 0.48 + 1.57461 -0.176<y ; 
Thus, the EOS proposed by Soave is of the form 

„ RT a(T) 
v-b v(v + b) 

In Equation (2.67) a can be expressed as 

R~T~ 

a = 0.42748 - 

P 


(2.68) 


(2.69) 


(2.70) 


1 1 


RT 

b. =0.08664— 

* P 

L 


(2.71) 



When calculating fugacity coefficients, it is convenient to write the equation in terms of 
the compressibility factor Z. In this form the Soave equation of state is 


Z- - Z 2 +(A- B - B 2 )z - AB = 0 

(2.72) 

where 


r i 

f — i 

<3 r ‘ 

II 

(2.73) 

RT 

(2.74) 

Z = — 

RT 

(2.75) 

The van der Waals one fluid mixing rules are used for calculating 

a and b of the mixture. 


Thus the fugacity coefficient of component i in a mixture is given by 

M , =— (Z-l)-ln(Z-fl) + — 

' b K ’ v ’ bRT 


~> 


In 


Z + 5 


(2.75) 


2.3.5 Advantages and Disadvantages of Equation of state approach 

In the Equation of State approach the definition of fugacity coefficient is used to 
describe both liquid and vapor phases through the application of single equation of state. 
The major advantage of this approach is its applicability over a wide range of pressures 
including critical and supercritical pressures. However, phase equilibrium calculations by 
these methods are rather lengthy except for a simple equation of state. Another problem 
associated with the use of the EOS technique is the insufficiency of information about the 
exact form of the EOS for mixtures, and the inaccuracy of the existing mixing rules for 
mixtures containing polar and hydrogen bonding components. 



CHAPTER S 

LITERATURE REVIEW 

The synthesis and design of separation processes requires a reliable knowledge of the 
phase behavior of the system to be separated. For the description of multicomponent non- 
electrolvte systems. g E -models or equations of state can be applied using binary data 
alone. Experimental data are however often missing. In such cases group contribution 
methods such as ASOG (Derr and deal. 1969:Kojima and Tochigi. 1979). UNIFAC 
(Fredenslund et al.. 1975.1977) and modified UNIFAC (Dortmund) (Weidlich and 
Gmehling. 1987:Gmehling et al.. 1993:Lohmann et al.. 1998:Larsen et al.. 1987) can be 
successfully applied. These methods are developed for the prediction of vapor-liquid 
equilibria (VLE). so that the required group interaction parameters were mainly fitted to 
experimental VLE data stored in the Dortmund Data Bank (Gmehlin.. 1985. 1991). The 
actual parameters of these methods were published by Hansen et al. (1991) and Tochigi 
etal. (1990). 

Koji Arita et al. ( 1991 ) measured the isothermal vapor-liquid Equilibria (VLE) data for 
the systems chlorofluoromethane (HCFC-22) and 1.1.1.2-tetrafluoroethane (HFC-134a) 
over the temperature range 273.16 to 323.16 K. The data obtained were correlated with 
use of the Redlich-Kwong-Soave (Soave. 1972) equation of state with one interaction 
parameter. The critical constants of HCFC22 and HFC134a were taken from the 
literatures (Japanese Association of Refrigeration. 1975: Wilson and basu. 1988). The 
acentric factors used in this study, which were calculated by using the generalized Frost- 
Kalkwarf-Thodos vapor-pressure correlation (Reid et al.. 1977). were 0.2222 for 
HCFC22 and .3259 for HFC134a. respectively. 

Hironobu Kubota et al. (1990) measured the isothermal VLE data for the system 
chlorodifluoromethane (HCFC22) + 1-chloro-l. 1 -difluoroethane (HCFC142b) over the 
temperature range 263 to 338 K and pressure 0.1 tO 2.7 MPa. Data at 263.15 and 293.15 
K measured by Nagahama et al.show maximum deviation of 2.2 and 3.3 mole percent 
with present experimental data respectively. To correlate the present phase equilibrium 
data, the Peng-Robinson equation (PR) of state with suitable binary interaction parameter 
kjj used. 
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VLE data for seven binary fluorocarbon mixtures were modeled using the Wong- 
Sandler (WS) model (1998). A variety of fluorocarbon azeotropes were chosen for testing 
the accuracy of the model. The azeotropes tested were CF 2 H 2 and CF 3 CH 3 . CF 2 H 2 and 
CF 3 CF 3 . CF 2 H 2 and CHF 2 CF 3 . CF 2 H 2 and CH 3 CH 2 CH 3 . CF 3 H and CF 3 CF 3 . CF 3 CF 3 and 
CO 2 . CHF 2 CF 3 and NH 3 . Using the WS mixing rule and non-random two liquid (NRTL) 
activity' coefficient model, a three parameter fit was calculated for each binary mix ture. 
The single adjustable parameter (kjj) for each binary pair is referred to as the WS mixing 
rule parameter (k. 12 ). Four excess free energy (or equivalently activity coefficient) 
models were evaluated for calculating G ex : UNIQAC. Van Laar. Wilson, and NRTL. In 
all cases the. NRTL model with WS mixing rule provided the best fit with VLE data. 

Yun Whan Kang et al. (1996) measured the isothermal VLE data for the three binary 
systems HFC-32 + HCFC-22. HFC-32 + CFC-12. and HFC-32 + HCC-40 at 10°C and 
correlated with the Wilson equation. The difluoromethane (HFC-32) + 

dichlorodifluoromethane (CFC-12) system forms a minimum boiling azeotrope, but the 
other do not. 

Min Young Jung et al. (2001) measured the isothermal VLE data of the binary system 
difluoromethane (HFC-32) + pentafluoroethane (HFC-125) over the temperature range 
from 268.15 K to 308.15 K and correlated with the Camahan-Starling-De Santis (CSD). 
Peng-Robinson (PR), and Redlich-Kwong-soave (RKS) equations of state. Higashi and 
Kobavashi and Nishiumi have published experimental data for the VLE for HFC-32 
+HFC-125. The binary interaction parameter for the CSD equation of state has stronger 
temperature dependence than those for the RKS and PR equations of state. The average 
deviation between the measured pressure and calculated results from the CSD equation of 
state is about 0.18% and that from the PR equation state is about 0.22%. The average 
deviation from RKS equation of state is about 0.58%. As a result. The CSD equation of 
state correlates the experimental data better than the PR and RKS equations of state. 

Eui-Yeop Chung et al. (1997) measured the isothermal VLE data of the binary mixture 
of HFC-32/134a over the temperature range from 263 K to 323 K and correlated with the 
Peng-Robinson and Redilch-Kwong-Soave equations of state. The relative pressure 
deviation for two cases was generally less than 0 . 2 % for the temperature range 
investigated, and the difference is not significant considering the uncertainty of pressure 
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investigated, and the difference is not significant considering the uncertainty of pressure 
measurements. Measured deviation of pressure and composition using Peng-Robinson 
equation of state correlates experimental data better than the Redlich-Kwong-Soave 
(RKS) equation of state. 

Chang Nyeon Kim et al. (2000) measured the isothermal VLE data of the binary system 
of difluoromethane (HFC-32) + 1.1.1-trifluoroethane (HFC-143a) over the temperature 
range from 263.15 K to 313.15 K and correlated with The Canahan-Starling-DE Santis 
(CSD). Peng-Robinson (PR), and Redilch-Kwong-Soave (RKS) equations of state using 
van der Waals mixing rules. It is evident that the binary interaction parameter for the 
CSD equation of state increases slightly as the temperature increases. On the other hand, 
the binary interaction parameters for both the PR and RKS equations of state decrease as 
the temperature increases. The binary interaction parameter for the RKS equation of state 
has larger temperature dependence than those for the CSD and PR equations of state. The 
average deviation between the measured pressure and the calculated results from the CSD 
equation of state is about 0.29% and that from the PR equation of state is aboutO.45%. 
The average deviation from the RKS equation of state is about 0.70%. As result, the CSD 
equation of state correlates the experimental data better than the PR and RKS equation of 
state. 

Jaewon Lee et al. (1996) measured the P. T. x. v data for HFC-134a +HCFC-124 and 
HCFC-124 + HCFC-142b and correlated the results with the Peng-Robinson equation of 
state using the van der Walls one fluid mixing rule. Neither HFC-134a +HCFC-124 nor 
HCFC-124 +HCFC-142b shows azeotropic behavior at the temperature studied. 

Sergio Bobbo et al. (1998) measured the VLE data for the systems 1. 1.1.2- 
tetrafluoroehane +isobutene (R134a + R600a) at 293.66 and 303.68 K. 1. 1.1.2- 
tetrafluoroethane + 1.1.1.3.3.3-hexafluoropropane (R134a + R236fa) at 283.62 and 
303.68 K and isobutene +1.1.1.3.3.3-haxafluoropropane (R600a + R236fa) and correlated 
with Redlich-Kwong-Soave equation of state with the Huron- Vidal mixing rules was 
used the relevant Margules equation. The R134a + R236fa system was almost ideal. On 
other hand both of the binaries R134a + R600a and R600a + R236fa composed of a 
hydrocarbon (here R600a) and the relatively highly fluorinated R134a or R236fa. 
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exhibited positive homogeneous azeotropes with qualitatively similar deviations from 
Raoult s law. 

C.N Kim et al. (1999) measured the isothermal VLE date of the system HFC-143a + 
HFC-134a over the temperature range from 263.15 to 313.15 K at 10 K intervals and 
correlated with the Camahan-Starling-De Santis (CSD) and Peng-Robinson (PR) 
equation of state. The Binary interaction parameters for both the PR and the CSD 
equations of state decrease as the temperature increases. The present data were in go'od 
agreement with the calculated results from the CSD equations of state, and the deviations 
were less than 1 .0% with the exception of two points. 

Michael Kleiber (1995) measured the VLE data for the systems R134a-R23. R134a- 
R114. R13B1-R23. R134a-dimethylether. R13B1-R13 and R142b-dimethvl ether over 
the temperature range 251 tO 298 K. and pressure up to 2MPa and correlated with 
UNIFAC equation. The consistency of the data has been shown by a maximum likelihood 
method. This paper also published optimized UNIFAC group assignment for refrigerants. 

Michael Kleiber (1994) and coworkers obtained the VLE for the binary mixtures 
propvlene-R12. propvlene-R22. propylene -R1 52a. propyl ene-Rl 3. propvlene-R23. 
propvlene-Rl 3 4a. propvlene-Rl 15. propvlene-R13Bl, propvlene-Rl 14. propvlene-Rll'6. 
R134a-R12. R134a-R152a. R134a-R116. R134a-R12. R134a-propane at temperatures 
between 25 1 and 298 K and pressure up to 2 MPa. He showed the consistency of the data 
by a maximum-likelihood method. He also correlated the data by generalized forms of 
various EOS such as Peng-Robinson. Peng-Robinson-Strvjeck-vera. Redlich-kwong- 
soave and Lee-Kesler-Plocker along with appropriate mixing rules. These EOS were 
fitted to the experimental data by optimizing the binary interaction parameters. 
Considering the fitting rules. The LKP equation seems to perform worse than the others, 
especially for azeotropic mixtures like propylene-R152a. propylene-Rl 15. propylene- 
R134a. R134a-propane R134a-R12.The PR. RKS and PRSV equations yield similar 
results in each case. The system R134a-Rl 16 was described significantly worse than the 
other mixtures. Only the PRSV equation with the Vanlaar type was able to describe it 
satisfactorily. This mixing rule is. of course, superior to the others, because it has two 
adjustable parameters instead of only one. but usually it gives only slightly better results 
than normal vdW mixing rule 



CHAPTER 4 


PREDICTION OF VLE DATA 


The purpose of phase-equilibrium thermodynamics is to predict the conditions (T. P. 
composition) prevailing when two or more phases are in equilibrium. The 
thermodynamic equations which determine the state of equilibrium between vapor and 
liquid phases are T 1 =T' (T is constant). P‘ = P' (P is constant), and f[ = . To fin d 
the conditions which satisfy these equations, it is necessary to have a method for 
evaluating the fiigacity ( f] ) of each component in both vapor and liquid phases. This can 
be done using the basic thermodynamic equation for calculating the fugacity coefficient 
( <j>\ or ft ). 

The fugacity coefficient can be obtained by using the relation 
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(4.2) 


The solution to the phase-equilibrium problem is provided completely by either one 
of these equations together with an EOS and the equations of phase-equilibrium. Eq. (4.1) 
is used when the EOS is given in volume explicit form ( V = V v (P.T . «,.«-,..)) whereas 

Eq. (4.2) is used when the EOS is in the pressure explicit form (P = P, (V.T . «,.«,...)) . 

Most of the EOS are pressure explicit type. To use an EOS. one must transform the EOS 
in molar form to a total form first. When one is to calculate the fugacity of component i 
in the mixture at a specified T. P. and composition, differentiate the EOS w.r.t n,. and 
substitute in Eq. (4.2) to get the result. 



There are basically four classes of VLE problems: 
BUBLP : Calculate (y, } andP. given {x,} . T. 


DEWP: Calculate {x,} and P. given [>•,}. T. 


BUBLP : Calculate {y,} and T. given {x,}. P. 
DEWT : Calculate {x,} and P. given [>•,}. P. 


The problems of engineering interest generally deal with BUBLP or DEWP. In order 
to correlate the experimental data, choice of the model to calculate the fugacitv 
coefficients and the determination of the best or the most representative parameters for 
the model are the most important steps. The models in this chapter are based on the PR 
EOS and Quartic EOS. each combined separately with vdW mixing rules and 
composition dependent two parameter mixing rules. 

4.1 UNIquac Functional group Activity Coefficients (UNIFAC) Method 
(group contribution method) 

To correlate and to predict the thermodynamic properties of a solution it was 
convenient to consider the component molecules as a collection of functional groups. 
The functional groups are structural units such as -CH3. -OH and others which when 
added form the constituent molecules. In the group contribution methods, a solution of 
components is treated as a solution of groups. The activity coefficients of the 
components are then determined by the properties of the functional groups rather than by 
those of the molecules. The basic idea behind this approach is that there are several 
thousands of chemical compounds that are of interest to a chemical engineer but the 
number of functional groups, which constitute these compounds, is very small. The 
activity coefficients in a large number of solutions can be calculated from the parameters 
characteristic of the functional groups. 
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In the UNIFAC method the activity coefficients consist of two parts, a combinatorial 
part (lny, c ) and a residual part (lny, R ). To evaluate Inyb only pure component data is 
required. 

c R 

lny, =lny, (combinatorial) + In r, ( residual ) (4-3) 

Inv," = In — + — q. In — + /, - — X-v./. (4.3) 

x, 2 x, , 

where /, = ^ (r, - q, ) - (r. - 1 ) (4'.4) 


The structural parameters r, and q, are calculated as the sum of product of the group 
volume and parameters Rk and Ok. 


q, 


(4,5) 


where v* is the number of groups of type A: in a molecule of component z. 
The residual part of the activity coefficient lny, R is given by 
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x„, = mole fraction of group m in the mixture 
a mn = group interaction parameter (in K). a mn # a nm 
u mn = measure of interaction energy between groups m and n 
Fk = group residual activity coefficient; and 

Fk" = residual activity coefficient of group k in a reference solution containing 
only molecules of type i. 


4.2 Modified UNIquac Functional group Activity Coefficients (Modified 
UNIFAC) Method (group contribution method) 


In the Modified UNIFAC method the activity coefficients consist of two parts, a 
combinatorial part (Inyf ) and a residual part (lny, R ). To evaluate Inyf only pure 
component data is required. 

l R * 

lny, = iny/ {combinatorial) + In;/, ( residual ) (4.10) 


The combinatorial part of the activity coefficient Iny, c is given by 
In;/ " =1-1’ + Ini’ -5q,\ 
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The parameter I' can be calculated by using the relative van der Waals volumes Rk 
of the different groups. 
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(4.15) 
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where v/" is the number of groups of type A: in a molecule of component i. 


(4.16) 


The residual part of the activity coefficient lny, R is given by 
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x,„ = mole fraction of group m in the mixture 
u mn = group interaction parameter (in K). a m „*a nm 
b m „ = group interaction parameter (in K). b mn *b nm 
c mn = group interaction parameter (in K). c m „^c nm 
u m n = measure of interaction energy between groups m and n 
fk = group residual activity coefficient: and 


(4.20) 
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rf = residual activity coefficient of group k in a reference solution containing 
only molecules of type i. 


4.3 Peng-Robinson Equation with van der Waals one fluid mixing rules 
4.3.1 Peng-Robinson Equation of State 

In 1976. Peng and Robinson introduced another equation of state, which extended 
some of the ideas presented by Soave. The repulsive part of the van der Waals model 
remained the same, but just as Soave had done, the temperature dependence of the 
attractive term was incorporated into the a term. The denominator of the attractive term 
was also changed to a slightly more complicated expression of volume and the size 
parameter. The Peng-Robinson Equation of State is given by 


p RT a(T) 

(4.21) 

v — b v(v + b) + b(v - b) 

RT 

where b = 0.07780 — - 

P 

L 

(4.22) 

R'-T 2 


a(T ) = 0.45724 

(4.23) 

a(T) = a(T )a 

(4.24) 

i i 

a 5 =\ + S(l-T}) 

(4.25) 


where reduced temperature 

(4.26) 

(4.27) 



S = f(co) = 0.37464 + 1 .54226<y - 0.26992<y : 
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and co is the Pitzer ( 1 995) acentric factor defined as 


ZW = — Iog| 0 P' [, _ =0 - -1.0 

(4.28) 

where P r ' s = ~ reduced pressure at T r =0.7 

(4.29) 

43.2 van der Waals one fluid Mixing Rules 


These are also called the classical mixing rules and are given by: 



(4.30) 

II 

-M' 

~M' 

(4.31) 

Where a is the energy parameter and b is the volume parameter. a„ 

and b„ are the cross 

terms i.e. the interaction terms when ( / * / ) whereas i and / represent the components(c). 

1 

sf 

ii 

cf 

(4.32) 

(b„ +b,. ) 

K= a -A,) 

V “ ) 

(4.33) 

Where k„ and I„ are the binary interaction parameters: /,, being zero most of the times. 

When Eq. (4.34) is applied to Eq. (4.33) with /„ = 0 . it reduces to a linear term. 

-cT 

w 

II 

(4.34) 


The above vdW one fluid mixing rules are used for van der Waals. Redlich-Kwong. 
Soave-Redlich-Kwong. and Peng-Robinson equations of state with a single binary 
interaction parameter k,,. (/,, = 0). 

The PR EOS expressed on a total basis, in pressure explicit form, is 

P = Jh*L ,4.35, 

I —n,b l (l +n,b) + n l b{l -n,b) 

where V=vn, and n, = total number of moles 
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Eq. (4.35) when differentiated on molar basis and substituted in Eq. (4.2) with the use of 
van der Waals one fluid mixing rules (Eqs. (4.30) through (4.34)) gives the following 
expression for the fugacity coefficient. 


M, = — (Z — l)-ln(Z — £)- 
b 


a 

h. _? IK 

In 

z+b[\+4i) 

2'JlbRT 

b ~ V ci 

Z + S(l-V2) 


(4.36) 


where 


RT 


(4.37) 


Eq. (4.36) is used to calculate the fugacity coefficient of either liquid or vapor 
phase, by using the appropriate compressibility factor and the composition. 

4.4 VLE Calculations 

The most important step of calculating the fugacity coefficients ( <j>[ or <j>J ) is carried 

out by using the Peng-Robinson equation. This requires Z ' . the liquid phase 
compressibility factor or Z' . the vapor phase compressibility factor. These are obtained 
by solving the EOS which is expressed in terms of the compressibility factor Z. 

For PR EOS. we can obtain the Z ! and Z' by solving the following cubic EOS 
which is expressed in terms of the compressibility factor Z. 


where 


Z ? +(R-l)Z : +(.4-3J5 ; -2B)Z + (B 2 + B : ' -AB) = 0 


A = 


aP 

R 2 T : 



z = 


PV 

nRT 


(4.38) 


(4.39) 

(4.40) 

(4.41) 
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4.4.1 Newton -Raphson Method. 


In general Neton -Raphson method expressed is 




/(*,-) 

fw 


(4.42) 


Above Eq. (4.38) expressed interms Newton Raphson method. 


f(Z) = Z : ' +(B-\)Z z +(A-3B 2 -2B)Z + (B z + B : ' -AB) 



/(Z,) 

f\Z,) 


(4.43) 


(4.44) 


above Eq. (4.44) solving to obtain one root of Zj orZ,' . To find the other two roots using 
Quadratic equation. 


4.4.2. Solution of Quadratic Equation. 

In general form of the Quadratic equation is 

ax 2 +bx + c = 0 (4.45) 

discriminant = b~ —4 ac (4.46) 

if disriminant is less than zero roots are imaginary, discriminant is equal zero roots are 
equal and discriminant greater than zero roots are real. Two roots are expressed is 


root\ = 


-b + \lb~ : -4 ac 


2a 


rool2 - 


-b-4b~ : -4 ac 


la 


(4.47) 


(4.48) 


Compressibility factor expressed in terms of Quadratic equation. 

Z : + ( Z{ - (1 - B) ) Z + ( Z, /: - Z\ ( 1 - B) + ( A - 35- - 25) ) = 0 (4.49) 

above Eq. (4.49) write in terms of general formula 

a=l. (4.50) 

b = [Z[ -(1-5)) (4.51) 

c = ( Z'; - z; ( 1 - B) + ( A - 3B 2 - 2B ) ) (4.52) 
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The values of a. b and c substituting Eq. (4.46) through (4.48) to obtain two roots of the 
liquid or vapor compressibility factors. 

Eq. (4.38) when solved gives three roots. There are two possibilities: either all the three 
roots are real or only one real root (remaining two being complex). When there are three 
real roots, the smallest root is the liquid phase compressibility factor and the largest root 
is the vapor phase compressibility factor, intermediate root being meaningless. If it is 
only one root, it is taken either as the liquid phase or vapor phase compressibility factor 
depending on the calculations. 


4.6 Data Reduction 


As mentioned earlier, there are two important steps i) model selection and ii) 
evaluation of the best parameters for the model. This needs reduction of the VLE data. In 
this work, only P-x-y data (given T and x,) is used. The average absolute deviations in 
pressure and composition are calculated from the relations. 




where. P exp . P ca i are experimental and calculated pressure respectively. y lexp . y ku/ are 

experimental and calculated vapor phase compositions of component 1. and N is the 
number of data points. 
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^ START ^ 

/ READ T .x,.Pf. 
/ v T-R k ~Qk-a„ m -a h 



Calculate combinatorial activity 7 coefficient 


In y; =ln^ + ^,ln^ + /,--^Xx ; /, : 

=z:(r, -q,)-{r, -l) • r , =1 q, =IvJ' ) 0 ,. 

2 a- * 



Fig 4.1. Algorithm for VLE calculation using UNIFAC Method. 


J J 









Fig 4.2. Algorithm for VLE calculation using Modified UNIFAC 
Method. 


34 










Fig 4.3. Algorithm for bubble pressure calculation using an 
equation of state 
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Subroutine PR 



Fig 4.4. Subroutine for Peng-Robinson equation of state. 
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CHAPTER 5 


RESULTS AND DISCUSSION 


5.1 Results 

The results obtained from the UNIFAC. modified UNIFAC and Peng-Robinson with van 
der Waals one fluid one mixing rule are presented in this section. The results obtained 
with these three models are compared to ascertain the best model. A total of 3 1 binary 
systems at different temperatures were studied. These systems are classified into 15 
different classes, such as HCFC/HFC. HCFC/HCFC. HFC/CFC. HFC/HCFC. HFC/HFC. 
HFC/FC. BrFC/CFC. HFC/HC. HFC/HOC. HCFC/HOC. HC/CFC. HC/HCFC. HC/HFC. 
HC/FC. HC/BrFC. Table (5.1) represents the source of the experimental data of all the 
systems and study all the systems different temperatures. The nomenclature, critical 
constants and acentric factors of the compounds (Reid et al. 1987) are listed in Table 
(5.2). Optimized UNIFAC group interaction, volume and surface area parameters of the 
refrigerants are listed in Table (5.3). Optimized Modified UNIFAC group interaction, 
volume and surface area parameters of the refrigerants are listed in Table (5.4). Only 
isothermal VLE (BUBLP) data is considered. 

The average absolute deviations in pressure and vapor-phase composition of the all 
systems studied are presented in Table (5.5). The pressure (P) and vapor-phase 
composition ( 17 ) of one system at one temperature are presented in Tables (5.6) - (5.27). 
The minimum percent deviations observed for pressure and vapor-composition are 0.49 
and 0.35 respectively and the corresponding maximum values are 30.80 and 14.52. 
Figures (5.1) - (5.11) represents the pressure (P). liquid composition (x). and vapor 
composition (v) of the each system at particular temperature. 


37 



5.2 Discussion 

In this work thirty-one systems of different classes at different temperatures have 
been studied using the UNIFAC. modified UNIFAC and PR EOS with vdW one fluid- 
mixing rule. The performances of the models are different for different classes of systems 
considered in the work. 

For HCFC / HFC class of system (R22-R134a). which is non-azeotropic in nature, 
the P-R EOS with vdW one fluid mixing rule Predicted the VLE data better than other 
two models. For HCFC / HCFC class of system (R22-R142b and R124-R142b). For R22- 
R142b system, the VLE data obtained by the P-R with vdW one fluid mixing rule better 
than other two models and R124-R142b system, the VLE data obtained by the UNIFAC 
model better than other two models. 

For HFC 1C FC class of systems (R32-R12 and R143a-R12). the VLE data obtained 
by the UNIFAC model better than other two models. For HFC/HCFC class of systems 
(R32-R22. R134a-Rl 24 and R134a-R142b). for R32-R22 system, the VLE data obtained 
by the UNIFAC model better than the other two models. R134a-R124 system, the VLE 
data obtained by the Modified UNIFAC model better than other two models and R134a- 
R142b system, the VLE data obtained by the P-R with vdW one fluid mixing rule better 
than other two models. 

For HFC / HFC class of systems (R32-R125. R32-R134a. R32-R143a. R134a- 
R152a. R134a- R236fa and R143a-R134a). the VLE data obtained by the P-R with vdW 
one fluid mixing rule better than other two models. For HFC /FC class of systems (R23- 
R1 16 and R134a-Rl 16). the VLE data obtained by UNIFAC model better than other two 
models. For BrFC / CFC class of system (R13B1-R13). the VLE obtained by the P-R 
with vdW one fluid mixing rule better results than other two models. 

For HFC ' HC class systems (R134a -R290 and R134 -R660a). the VLE data 
obtained by the Modified UNIFAC model better results than other two models. For 
HFC/HOC class systems (R134a-Dimethyl ether), the VLE data obtained by the 
Modified UNIFAC model better than other two models. . For HCFC/HOC class systems 
(R142b-Dimethyl ether), the VLE data obtained by the Modified UNIFAC model better 
than other two models. 
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For HC / CFC class of systems (Propylene-R12. Propvlene-R13 and Propylene- 
R115). the VLE data obtained by the UNIFAC model better than other two models. For 
HC / HCFC class of systems (Propvlene-R22 and Propvlene-R142b). the VLE data 
obtained by the UNIFAC model better than other models. For HC /HFC class of systems 
(R600a-R236fa. Propvlene-R134a and Propvlene-R152a). the VLE data obtained by the 
UNIFAC model better than other two models. For HC / FC class of system (Propylene- 
R116). the VLE data obtained by the UNIFAC model better than other two models. For 
HC / BrFC class of system (Propvlene-R13Bl). the VLE data obtained by the Modified 
UNIFAC model better than other two models. 

This study indicates that well known UNIFAC model displays better abilities in 
general as compared to the Modified UNIFAC model and Peng-Robinson EOS with van 
der Waal one fluid-mixing rule. But for HCFC/HFC. HFC / HFC class of systems P-R 
with vdW one fluid-mixing rule gave better results and HFC/HC. HFC/HOC. 
HCFC/HOC .class of systems Modified UNIFAC method gave better results. 
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Table 5.1. Source of Experimental data for the systems. 


S.No 

System 

Data Source 

1. 

R22 - R134a 

Koji Arita et al. (1991) 

9 

R22 - R142b 

Hironobu Kubota et al. ( 1 990) 

*> 

J. 

R23-R116 

M.B. Shiflett et al. (1998) 

4. 

R32-R12 

Yun Whan Kang et al. (1996) 

5. 

R32-R22 

Yun Whan Kang et al. ( 1 996) 

6. 

R32-R125 

Min Young Jung et al. (2001) 

7. 

R32 - R134a 

Eui-Yeop Chung et al. (1997) 

8. 

R32 - R143a 

Chang Nyeon Kim (2000) 

9. 

R13B1 - R13 

Michael Kleiber (1995) 

10. 

R124 — R142b 

Jaewon Lee et al. (1996) 

11. 

R134a - R12 

Michael Kleiber ( 1 994) 

12. 

R134a - R1 16 

Michael Kleiber (1994) 

13. 

R134a - R124 

Jaewon Lee et al. (1996) 

14. 

R134a-R142b 

Michael Kleiber (1994) 

15. 

R134a-R152a 

Michael Kleiber (1994) 

16. 

R134a-R236fa 

Sergio Bobbo et al. (1998) 

17. 

R134a-R290 

Michael Kleiber (1994) 

18. 

R134a-R600a 

Sergio Bobbo et al. (1998) 

19. 

R134a - Dimethyl ether 

Michael Kleiber (1995) 

20. 

R142b - Dimethyl ether 

Michael Kleiber (1995) 

21. 

R143a- R134a 

C.N. Kim et al. (2000) 

22. 

R600a - R236fa 

Sergio Bobbo et al. (1998) 

23. 

Propylene -R1 2 

Michael Kleiber (1994) 

24. 

Propylene - R13 

Michael Kleiber ( 1 994) 
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(Continued) 


25. 

Propylene - R22 

Michael Kleiber ( 1 994) 

26. 

Propylene -R13B1 

Michael Kleiber (1994) 

27. 

Propylene - R1 1 5 

Michael Kleiber (1994) 

28. 

Propylene - R1 16 

Michael Kleiber (1994) 

29. 

Propylene - R134a 

Michael Kleiber (1994) 

30. 

Propylene - R142b 

Michael Kleiber (1994) 

31. 

Propylene - R1 52a 

Michael Kleiber (1994) 
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Table 5.2. Nomenclature, Critical properties and acentric factors of the refrigerants. 
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Table 5.3. Optimized UNIFAC group assignment for 
refrigerants 


Main 

group 

Subgroup 

No. 

Rt 

Qk • 

1 

ch 2 

ch 3 

1 

0.9011 

0.848 



ch 2 

2 

0.6744 

0.540 * 



CH 

j 

0.4469 

0.228 \ 



C 

4 

0.2195 

0.000 

40 

cf 2 

cf 3 

74 

1.4060 

1.380 • 



cf 2 

75 

1.0105 

0.920 



CF 

76 

0.6150 

0.460 

45 

cf 4 

cf 4 

86 

1.8016 

1.840 

51 

CHF 

ch 3 f 

109 

1.2966 

1.308 



ch 2 f 

110 



1.0699 

1.000 



CHF 

Ill 

0.8420 

0.688 

52 

chf 2 

ch 2 f 2 

112 

1.4654 

1.460 ; 



chf 2 

113 

1.2380 

1.232 



ch 3 - chf 2 

114 

2.1391 

2.080 



(C F 3 -)CH2F 

115 

1.0699 

1.000 

53 

chf 3 

chf 3 

116 

1.6335 

1.608 



(CH 2 -) cf 3 

117 

1.4060 

1.380 • 

54 

CC1F 

cci 3 f 

118 

3.0356 

2.644 



cci 2 f 

119 

2.2446 

1.920 

55 

CC1F 2 

cci 2 f 2 

120 

2.5926 

2.368 



ccif 2 

121 

1.8016 

1.644 



CBrF 3 

122 

2.4028 

2.232 

56 

CCIFs 

CC1 f 2 

123 

2.1971 

2.104 



(CF 3 -)CC1F 2 

124 

1.8016 

1.644 



CClBrF 2 

125 

2.4028 

2.232 
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(Continued) 


57 

CHC1F 

CHC1F: 

126 

2.0290 

1.872 

(CH 3 -) CC1F: 

127 

1.8016 

1.644 

CHC1F 

128 

1.6335 

1.412 

58 ^ 

CHCbF 

CHCbF 

129 

2.4562 

2.144 

CH 3 - CCFF 

130 

3.1457 

2.768 

CHCb- CF 3 

131 

3.6624 

3.378 

59 

CF 3 CI 3 

CFCb- CF 2 CI 

132 

4.0462 

3.564 
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Table 5.4. Optimized Modified UNIFAC group assignment for 
refrigerants 


Main 

group 

Subgroup 

No. 

Rk 

Qk 

1 

ch 2 

ch 3 

1 

0.6325 

1.0608.. 



ch 2 

2 

0.6325 

0.7081 



CH 

-* 

j 

0.6325 

0.3554 



C 

4 

0.6325 

0.0000 

40 

cf 2 

cf 3 

74 

1.7874 

1.7975 



cf 2 

75 

1.0843 

1.1801 



CF 

76 

0.3806 

0.5841 

70 

cf 4 

cf 4 

148 

1.5375 

1.9920 

71 

CHF 

ch 3 f 

149 

1.2962 

1.3087 



ch 2 f 

150 

1.0685 

1.0007 



CHF 

151 

0.9386 

0.9700 

72 

chf 2 

CH 2 F 2 

152 

1.7529 

1.8147 



chf 2 

153 

1.1396 

1.4390 



ch 3 - chf 2 

154 

1.9820 

2.3811 



(C F 3 -)CH2F 

155 

0.7108 

1.1005 

73 

chf 3 

chf 3 

156 

1.7272 

1.9637 



(CH 2 -) cf 3 

157 

1.4246 

1.3487 

74 

CC1F 

cci 3 f 

158 

2.6679 

2.7843 



cci 2 f 

159 

2.0418 

2.0035 

75 

CC1F 2 

cci 2 f 2 

160 

2.3591 

3.0749 ' 



ccif 2 

161 

1.3574 

1.7688 • 



CBrF 3 

162 

2.2767 

2.8122 - 

76 

CCIFs 

CC1 f 2 

163 

2.1347 

2.0935 ■ 



(CF 3 -) CC1 f 2 

164 

1.0782 

1.0584 



CClBrF 2 

165 

1.9280 

1.9139 

77 

CHC1F 

CHC1F 2 

166 

1.8552 

2.0267 



(CH 3 -) ccif 2 

167 

1.6803 

1.5016 

. 
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(Continued) 




CHC1F 

168 

1.4422 

1.6403 

78 

CHC1 2 F 

CHCFF 

169 

2.2504 

2.9420 

CH 3 - CCbF 

170 

3.1638 

2.8374 

CHCF- cf 3 

171 

2.2703 

3.5910 

79 

cf 3 c1 3 

CFCI 2 -CF 2 CI 

172 

3.3963 

4.3036 
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Table 5.5. Average deviation in vapor - phase composition, M and pressure, j for systems 











0.0131 0.0054 0.0172 0.0189 0.0535 



(Continued) 





Table 5.6: VLE Data for R22 (1) - R134a (2) 
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Table 5.7: VLE Data for R22 (1) - R142b (2) 
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Table 5.8: VLE Data for R124 (1) - R142b (2) 
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Table 5.9: VLE Data for R134a (1) - R12 (2) 
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Table 5.10: VLE Data for R134a (1) - R124 (2) 
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Table 5.11: VLE Data for R134a (1) - R142b (2) 
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Table 5.12: VLE Data for R32 (1) - R125 (2) 
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Table 5.13: VLE Data for R134a (1) - R152a (2) 




Table 5.14: VLE Data for R143a (1) - R134a (2) 
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Table 5.17: VLE Data for R134a (1)-R600a (2) 
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Table 5.18: VLE Data for R134a (1) - Dimethyl ether (2) 






0.0189 0.0118 0.0249 



Table 5.20: VLE Data for Propylene (1) - R12 (2) 




Table 5.21 : VLE Data for Propylene (1) - R1 15 (2) 







Table 5.23: VLE Data for Propylene (1) - R142b (2) 




Table 5.23: VLE Data for Propylene (1) - R134a (2) 




Table 5.25: VLE Data for Propylene (1) - R152a (2) 

Temperature: 255.0 K 





Table 5.26: VLE Data for Propylene (1) - R1 16 (2) 




o 

o 


CM 


• * 

u 

3 

3 

U 

OJ 

CU 


s 

H 



CN 



Table 5.27: VLE Data for Propylene (1) - R13B1 (2) 

Temperature: 251.0 K 
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Figure 5.1. P-x-y diagram for R22 (1) / R134a (2) system at T-273.15 K 
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igure 5.2. P-x-y diagram for R22 (1) / R142b (2) system at T=263.15 K 
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Figure 5.3. P-x-y diagram for R134a (1) / R142b (2) system at T=268.0 K 
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Figure 5.4. P-x-y diagram for R13B1 (1) / R13 (2) system at T=255.0 K 
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Figure 5.5. P-x-y diagram for R32 (1) / R134a (2) system at T=263.15 K 
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Figure 5.7. P-x-y diagram for R142b (1) / Dimethyl ether (2) system at T=268.0 K 





Figure 5.8. P-x-y diagram for Propylene (1) / R12 (2) system at T=258.0 K 





Figure 5.9. P-x-y diagram for propylene (1) / R22 (2) system at T=258.0 K 
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Figure 5.10. P-x-y diagram for Propylene (1) / R152a (2) system at T=255.0 K 





Mole fraction (x, y) of R13B1 

Figure 5.11. P-x-y diagram for Propylene (1) / R13B1 (2) system at T=251.0 




CHAPTER 6 


CONCLUSIONS 

The activity coefficient is directly related to molar Gibbs energy. A total of thirty- 
one binary mixtures of different classes at different temperatures, have been studied using 
UNIFAC. modified UNI FAC and Peng-Robinson EOS with van der Waal one fluid 
mixing rule models and predicted results are compared with experimental data. It was 
found that the UNIFAC model is able to predict the VLE data much better than the other 
two models. 

For HCFC'HFC. HFC/HFC. BRFC/CFC class of systems and R22-R142b. R134a- 
R142b systems, the VLE obtained by the P-R with vdW one fluid mixing rule gave better 
results than other two models. For HFC/HC. HFC/HOC. HCFC/HOC. HC/BrFC class of 
systems and Rl34a-R124 system, the VLE data obtained by the Modified UNIFAC 
method gave better results than other two models. 

Further studies can be carried out for the prediction of Vapor-Liquid Equilibria data 
of a large number of systems using both EOS and Activity coefficient method with 
different mixing rules. 
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